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ABSTRACT 

We investigate the internal structure of clusters of galaxies in high-resolution N-body 
simulations of 4 different cosmologies. There is a higher proportion of disordered clus- 
ters in critical-density than in low-density universes, although the structure of relaxed 
clusters is very similar in each. Crude measures of substructure, such as the shift in 
the position of the centre-of-mass as the density threshold is varied, can distinguish 
the two in a sample of just 20 or so clusters; it is harder to differentiate between 
clusters in open and flat models with the same density parameter. Most clusters are 
in a quasi-steady state within the virial radius and are well-described by the density 
profile of Navarro, Frenk & White (1995). 

Key words: Galaxies: clustering — Cosmology: miscellaneous 
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1 INTRODUCTION 

The Virgo Consortium is an international collaboration 
whose aim is to investigate the large-scale structure of the 
universe via high-resolution N-body, hydrodynamical simu- 
lations. The ultimate goal is to relate the observed galaxy 
distribution to that of the matter fluctuations in the early 
universe and in particular to determine the biases between 
the two, both in position and velocity, on various scales. 

As a first step, we have carried out pure N-body simu- 
lations using a large number of particles (256^ ~ 17 million) 
and four different Cold Dark Matter (CDM) cosmologies. 
Other papers will report on the clustering evolution of dark 
matter and dark matter halos, the topology of large-scale 
structure, strong and weak gravitational lensing properties, 
etc. . 

In this paper we concentrate on the structure of galaxy 
clusters in our simulations. Because of our high-resolution, 
we can trace both the large-scale topology and the internal 
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structure of the clusters simultaneously, finding a smaller 
variation in cluster properties between difi'erent cosmologies 
than are found in simulations of isolated clusters. 

The main reason for studying the structure of galaxy 
clusters is to try to discriminate between difi'erent cosmo- 
logical models. In critical-density universes clustering con- 
tinues to grow to the present day, whereas in low-density 
universes it begins to decline after a redshift z ~ — 1. 
This means that clusters in low-density universes are ex- 
pected to be dynamically more relaxed and to have less sub- 
structure, steeper density profiles and rounder iso-density 
contours. Early results from simulations of isolated clusters 
showed strong variations of all three of the above properties 
with the cosmological density parameter, fio (e.g. Mohr et 
al. 1995, hereafter M95). We find a much weaker cosmo- 
logical dependence. Although we do indeed find a higher 
proportion of disordered clusters in our critical-density sim- 
ulations, the properties of relaxed clusters are very similar 
in all the cosmologies we examined. 

M95 found much steeper density profiles for clusters in 
low-density cosmologies than in critical-density ones. Cen 
(1994) arrived at a similar conclusion, but the evidence is 
rather weak and may instead refiect a dependence on mass: 
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there are many more massive clusters in his critical-density 
simulation than in his open ones. Subsequent studies (e.g. 
Jing et al. 1995, Crone et al. 1997) have shown that, when 
measured at radii enclosing the same mean overdensity, the 
slope of the density profiles of all clusters are fairly similar 
and, although they tend to be slightly steeper in low-density 
universes (i.e. the clusters are more isolated), this difference 
is unlikely to be detectable. 

M95 also obtained very few elongated clusters in their 
N-body, hydrodynamic simulations of low-density universes 
but other studies find most clusters to be flattened in all 
cosmologies. Splinter et al. (1997) found no change in shape 
with Qo- Jing et al. (1995) found a slight, but significant 
tendency for rounder clusters in low-density simulations, and 
the present work confirms this result. 

The velocity dispersion profiles of clusters are another 
potential discriminant of the density parameter. Jing & 
Borner (1995) find that the velocity dispersion declines 
with radius significantly more rapidly in low-density than 
in critical-density universes (however the effect is hard to 
measure even with 100 redshifts per cluster). By contrast, 
we find that dynamically relaxed clusters in all our simula- 
tions are well-fit by the model of Navarro, Frenk & White 
(1995), in agreement with Huss, Jain & Steinmetz (1997). 
Moreover, we do not find a strong cosmological dependence 
for the anisotropy parameter as reported in Crone, Evrard 
& Richstone (1994). 

Substructure (Richstone, Loeb & Turner 1992; Bartel- 
mann, Ehlers & Schneider 1993; Lacey & Cole 1993) seems 
to be the best indicator of ilo- Other measures merely re- 
fiect the different abundance of relaxed clusters whose prop- 
erties are relatively independent of the cosmological setting. 
Pinkney et al. (1996) applied a wide variety of tests to simu- 
lations of merging and isolated clusters and concluded that 
there was no single statistic which was guaranteed to de- 
tect it: instead they recommend a battery of a dozen tests, 
depending upon the data available and the type of substruc- 
ture it is desired to detect. 

Crone, Evrard & Richstone (1996) found that a good 
test for discriminating between different cosmologies is the 
shift in the centre of mass of the matter contained within 
an iso-density contour as the co ntou r level is varied and 
we concur with this (see Section 3.3). Other powerful in- 
dicators are moments of the squared-density distribution 
(Dutta 1995) or of the two-dimensional potential (Buote & 
Tsai 1995) or of the surface density (Wilson, Cole & Frenk 
1996). One promising new method for detecting substruc- 
ture is the 'hierarchical clustering method' which looks for 
dynamically bound associations of galaxies (Serna & Gerbal 
1996, Gurzadyan & Mazure 1997): however, this is much 
more complicated to apply than the above tests and so we 
do not consider it here. 

The above results are somewhat contradictory and con- 
fusing for several reasons: (i) variable mass resolution, (ii) 
different box-sizes: some clusters are isolated and some em- 
bedded within large-scale structures, (iii) different power 
spectra and methods of generating initial conditions, (iv) 
different mixtures of gas and dark matter, etc.. The data 
presented here provide a large sample of clusters covering 
about a decade in mass, drawn from large-scale simulations 
yet with a minimum resolution of more than 1000 particles 
within the virial radius. We feel that our cluster sample is 



Table 1. Parameters of the 4 simulations, as described in the 
text. 



Label 


^0 


Ao 


r 


h 




N 


m/{h-^ Mq) 


OCDM 


0.3 


0.0 


0.21 


0.7 


0.85 


200^ 


1.4 X 10" 


ACDM 


0.3 


0.7 


0.21 


0.7 


1.30 


256^ 


6.8 X lO^f 


tCDM 


1.0 


0.0 


0.21 


0.5 


0.68 


256^ 


2.3 X 10" 


SCDM 


1.0 


0.0 


0.50 


0.5 


0.61 


256^ 


2.3 X 10" 



much more homogeneous and representative of the real clus- 
ter population than those of previous studies. 

The outline of the paper is as follows. The simulations 
and numerical method are described in Section ^, as is the 
method of cluster identification. The properties of the clus- 
ters and various statistical measures that describe them are 
presented in Section ^ and the prospects for discriminating 
between various cosmologies are assessed in Section ^ 



2 METHODOLOGY 
2.1 The simulations 

We report on the results of 4 simulations in this paper. Their 
properties are chosen to be as similar as possible, while span- 
ning a range of cosmological parameters as listed in Table |l| 
The box-size is fixed at 239.5/i~^ Mpc in each case, where 
the Hubble parameter, h = ffo/lOOkms"^ Mpc~^ , is cho- 
sen so as to give a sensible age for the Universe (ranging 
from 11.5 Gyr for OCDM to 13.8 Gyr for ACDM). The power 
spectrum is taken to be that of cold dark matter (CDM), 
where for three of the runs the shape parameter, F, is set 
equal to the physically-motivated combination ^loh. In the 
fourth run, ACDM, we take F = 0.21 so as to give the same 
spectral shape as in the low-density models; this set of pa- 
rameters could occur, for example, in a decaying-neutrino 
model (Efstathiou, Bond & White 1992) 

It was our intention to choose the normalisation of the 
power spectrum, erg, so as to give approximately the same 
number of Abell clusters in each run (see, for example. 
White, Efstathiou & Frenk, 1993; Viana & Liddle, 1996; 
Eke, Cole & Frenk, 1996). In this we were only partially 
successful with a range of a factor of four in abundance; in 
fact, the normalisation is too large in all runs except OCDM. 
We have since repeated the other three runs with the correct 
normalisation but we choose to work with the original sim- 
ulations in this paper as this gives a larger cluster sample: 
the properties of the individual clusters (as contrasted with 
their large-scale distribution) are almost independent of ag. 

In run OCDM, the number of particles was only half 
that in the other simulations. Note, however, that the par- 
ticle mass is intermediate between that of the other low- 
density run, ACDM, and the critical-density runs. 

The gravitational softening in each run was different 
(because of the different mass-resolution) but was never 
larger than 30/i~^ kpc. Here we will limit our discussion to 
scales greater than 60h~^ kpc. 



2.2 The numerical method 

The initial conditions for the simulations were created by 
perturbing particles from a glass-like initial state (White 
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1996; see also Couchman, Thomas & Pearce 1995). The 
waves were drawn with random phases and amplitudes from 
a Gaussian power-spectrum (Bond & Efstathiou 1984) using 
the prescription of Efstathiou et al. (1985). 

The initial state was evolved forward from z ~ 50 
using the Hydra code of Couchman, Thomas & Pearce 
(1995), the parallelisation of which is discussed in Pearce 
& Couchman (1997). The timestep was set equal to At = 
min(0.25dta, 0.5dt„), where dta = mini(s/ai)^^^, dtv = 
mhii(s/vi), s is the softening, ai the acceleration, Vi the 
speed, and the subscript i runs over all particles. The num- 
ber of timesteps varied from 956 for the rCDM run to just 
under 1305 in run ACDM that has a smaller particle mass 
and softening. Energy non-conservation in each run, as mea- 
sured by the divergence in the Layzer-Irvine integral, was 
approximately 1 percent. 

2.3 Cluster identification 

Observationally, clusters of galaxies are rather poorly- 
defined objects. The traditional optical selection criterion 
depends upon the surface density of galaxies, but this is 
prone to projection effects. In this paper we adopt a more 
theoretical viewpoint, defining clusters as sets of particles 
which exceed a given overdensity threshold. 

The usual way to impose a density cut is to group 
together particles which lie within a certain distance, or 
'linking length', / (the 'friends-of-friends' approach: Davis 
et al. 1985). If this is set equal to h times the mean inter- 
particle separation, ie I = bn~^^^ , then the associated over- 
density is approximately 2&~^. In this paper we adopt a 
different approach whereby the density at the location of 
each dark matter particle is calculated by smoothing over 32 
neighbours using a smoothed particle hydrodynamics (SPH) 
algorithm (see e.g. Hernquist & Katz 1989). This has the 
advantage of giving a continuous range of densities whose 
values are much less subject to Poisson fluctuation than the 
nearest-neighbour approach. 

To identify clusters we first select those particles whose 
density exceeds some threshold, namely 180 times the criti- 
cal density. In order to treat the various cosmologies in the 
same manner, we use only a fraction r2o(256^/A*') of the par- 
ticles, chosen at random, at this stage (subsequent analysis 
uses all the particles). This ensures that each selected par- 
ticle corresponds to the same mass, 2.26 x 10^^h~^ M©. We 
then create a minimal spanning tree of the data which can 
be truncated at a range of linking- lengths, such as that given 
above, to divide the particles into groups. In practice, any 
linking length between and a few times this gives 

the same division, except where two groups happen to lie 
very close together. 

We define clusters in terms of an overdensity relative 
to the critical density because we wish to compare similar 
objects in each simulation. Eke, Cole & Frenk (1996) show 
that collapsed structures, in the spherical top-hat approx- 
imation, will have lower overdensities (relative to critical) 
in universes with Qq < 1. As the true value of the density 
parameter is unknown, however, we have no way of allowing 
for this in the observations. 

In principle we could search for groups down to rela- 
tively small masses. In this paper, however, we are inter- 
ested only in rich clusters and so we choose a mass-cut of 



Table 2. Cluster data: Afdus — the number of clusters with mass 
exceeding 2 X 10^'^ Mq, a — the mean 1-D velocity dispersion 
of a 2 X 10"h-i M0 cluster. 

Label A^clus cr/ kms~^ 

OCDM 90 720 

ACDM 208 730 

tCDM 377 690 

SCDM 289 680 



T 




O I , , , , 1 , , , , 1 , , L 

5x10^'^ 10^^ 

MQSs(6>180)/10^'^h"^M 

Figure 1. The mass of clusters in the SCDM cosmology defined 
by two different techniques: abscissa, mass within an overden- 
sity contour of 180 relative to the critical density; ordinate, mass 
within an Abell radius. 



2 X 10^'* /i^^ M0. The number of clusters in each simulation 
using this definition is listed in Table |^. Also shown is the 
mean 1-D velocity dispersion of 2 x 10^*/i~^ M0 clusters: 
these are similar in the different cosmologies which suggests 
that our selection criterion is a sensible one. The correspond- 
ing kinetic temperature is approximately 3 keV. The velocity 
dispersions are slightly higher in the low-density universes 
reflecting the fact that clusters in these cosmologies tend 
to form at higher redshift and hence be smaller and have a 
higher velocity dispersion for a given mass. 

An alternative way to deflne clusters is in terms of the 
mass contained within a sphere of fixed radius. We find a 
close correlation between the mass of clusters defined by the 
overdensity technique described above and the total mass 
within an Abell radius centred on the cluster centroid, which 
suggests that there is little difference between the two selec- 
tion methods. Figure |l| shows the correspondence for SCDM; 
the other cosmologies show similar behaviour. 

Throughout this paper we use the term 'the mass of 
the cluster' and the symbol miso as synonyms for the mass 
contained within an overdensity contour of 180 relative to 
the critical density. 



3 RESULTS 
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Figure 3. Values of the Mext statistic for clusters OCDM and 
ACDM, combined. 



Figure 3 - continued Values of the Mext statistic for clusters 
rCDM and SCDM, combined. 



3.1 Multiplicity function 

Theoretically, one of the most basic differences between clus- 
ters in critical and sub-critical density universes is that the 
latter should be more relaxed, having formed at a higher 
redshift. This should manifest itself in various ways, one be- 
ing that these clusters should be more isolated and have 
less substructure. Visually there is little difference between 
clusters in the different cosmologies: Figure ^ shows some 
plots of clusters of mass 3-4x 10"/i"^ Mq in the OCDM and 
rCDM runs. Consequently, we define a statistic to quantify 
the degree of isolation of a cluster. 

Firstly we catalogue all groups of particles with over- 
density greater than 180 times the critical density and which 
lie within 2 Abell radii of the cluster centre (this catalogue 
is obtained as a by-product of the method we use to identify 
clusters). Then we define a multiplicity statistic. 



Moxt = — 



(1) 



where M is the total mass of all the groups and m is the 
mass of the cluster. 

The majority of clusters are relatively isolated, with 
A^oxt < 1-2. However, there is a small fraction of clusters 
with larger values of Mext and this fraction increases with 
decreasing mass. The two low-density simulations have iden- 
tical distributions, as measured by a Kolmogorov-Smirnov 
test, as do the two critical-density universes. Moxt is plotted 
against mass for each pair of simulations in Figure ^ and 
a comparison of the two probability density distributions is 
shown in Figure ^. The critical-density simulations show a 
significantly higher fraction of clusters (approximately one 
quarter) which have neighbouring structures with a mass 
exceeding 20 percent of that of the cluster itself. Note, how- 
ever, that there are few binary clusters (for which at least 
one of the pair would have Mext > 2). We have checked that 
the differences in Mext are not due to the different mass dis- 
tributions of the two samples but are in fact present at both 
low and high masses. 
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Figure 4. Comparison of the probability distributions of the 
Afoxt statistic for clusters in low-density (solid line) and critical- 
density (dotted line) simulations. All the clusters with Mext> 2 
are placed in the right-hand bin. 



3.2 Substructure 

As well as being more isolated, clusters in low-density uni- 
verses are expected to show less substructure. We measure 
this by another multiplicity statistic. For each cluster we 
start with the merger tree used in its definition and then 
gradually decrease the linking length, I. This causes the clus- 
ter to break up into subclumps of higher-and higher over- 
density. At each stage we define 



M„ 



mi + 7712 + ma 
mi 



(2) 



where mi > m2 > ma are the masses of the three largest 
groups. The maximum value of Mint as I is lowered by a 
factor of 4.64 (corresponding to a density increase of a factor 
of 100) is used as a measure of the number of subclumps (up 
to a maximum of 3) in the catalogue. We experimented with 
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Figure 2. Example clusters from the OCDM simulation. The circles represent an Abell radius. 
© 0000 RAS, MNRAS 000, 000-000 




Figure 2 - 



continued Example clusters from the tCDM simulation. The circles represent an Abell radius. 
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Figure 5. Comparison of tlie probability distributions of the Afint 
statistic for clusters in low-density (solid line) and critical-density 
(dotted line) simulations. 



Figure 6. Comparison of the probability distributions of the C2 
statistic for clusters in low-density (solid line) and critical-density 
(dotted line) simulations. The second bin for Qq = 0.3 extends 
to 6.0; all values of C2 greater than 1.0 are garnered into the 
right-hand bin. 



slightly different definitions for the multiplicity statistic but 
all gave similar results. 

Once again the two Qq ~ 0.3 and the two flo ~ 1.0 
universes give similar distributions, this time with no mass- 
dependence. The probability distributions of the two pairs 
are compared in Figure The difference between the two is 
significant with the latter showing more substructure. The 
mean values of Mint are 1.28 and 1.53, respectively, and clus- 
ters in the critical-density models have twice the frequency 
of multiple cores (approximately 40 percent with Afint>1.5 
and 20 percent with Mint> 2). 

In combination, the Mext and Mint statistics show that 
there are very real differences in the degree of substructure 
between clusters in low-density and critical-density cosmolo- 
gies. The question remains as to whether these differences 
would be detectable in practice. Crone, Evrard & Richstone 
(1996) suggested that the shift in the position of the centre- 
of-mass of matter contained within an iso-density contour, 
as the level of the contour is altered, would be a good dis- 



radius, rigo (defined in Section 3.3), is shown in Figure n for 



overdensities in the range 180-18 000. With a sample of just 
20 clusters, these two distributions could be distinguished 
with 5 percent confidence. A very similar result is obtained 
by comparing the position of the centre-of-mass to that of 
the density maximum. Of course, the density contours of 
the dark-matter distribution are not directly measurable, 
but will be related to those of the X-ray emitting gas in a 
way which can be quantified by future simulations. 

If one is restricted to observations of galaxies only, then 
the best measures of substructure may well be the hier- 
archical clustering methods of Serna & Gerbal (1996) and 
Gurzadyan & Mazure (1997). The latter find that the exis- 
tence of 2-3 such subclumps is typical in a sample of clusters 
drawn from the ESO Nearby Abell Cluster Survey. Unfor- 
tunately their method is complicated to apply and will be 
deferred to a future paper. 



3.3 Morphology 

Continuing our theme of looking for evidence of relaxation 
in clusters in the low-density cosmologies, we next examine 
their shapes. M95 found enormous differences between the 
axial ratios in low-density and critical-density universes in 
their models with the former being much closer to unity. We 
define the semi-axes, ai > 02 > 03, of each cluster in terms 
of the best-fitting ellipsoid, i.e. Ui = (SAi)"'"" where the \i 
are the eigenvalues of the inertia tensor. 



(3) 



and the sum extends over all m particles in the cluster (the 
normalisation is chosen so that for a uniform sphere the 
semi- axes are equal to the radius). The mean radius at which 
the overdensity equals 180 is rigo = (010203)^^^^. 

Yet again the distribution of axial ratios in the two low- 
density simulations and in the two critical-density simula- 
tions are indistinguishable. The latter are significantly more 
elongated (as measured by a KS test) but only by a small 
amount. Figure ^ shows a histogram of the minor/major axis 
ratio for the two sets of runs. The mean value of the ratio 
is 0.50 in the low-density runs and 0.43 in the high-density 
runs. Thus there is a slight tendency for the clusters in low- 
density universes to be rounder than those in critical-density 
universes, that is consistent with them being slightly more 
isolated. There is also a weaker trend for high-mass clusters 
to be rounder than low-mass ones. 

The three-dimensional shapes of the clusters in all cos- 
mologies show a wide range of triaxialities, ranging from 
highly-prolate to highly-oblate. Figure ^ shows the distribu- 
tion in ellipticity-prolateness space for SCDM (all are similar 
in appearance). The ellipticity is defined as 



A3 ^ — Aj^ 



(4) 
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0.2 0.4 0.6 0.8 

axial ratio (min/max) 

Figure 7. Values of the minor/major axis ratio. = 0.3 refers 
to the combined sample ACDM and OCDM; VIq = 1.0 to SCDM 
and tCDM. 




0.5 



Figure 8. Distribution in the ellipticity-prolateness plane of clus- 
ters in the SCDM cosmology. 



(5) 



and the prolateness as 

Given that the contours of X-ray emission will be 
rounder than the iso-density contours, it is extremely un- 
likely that observations will be able to distinguish between 
low-density and critical-density universes based simply on 
the degree of flattening of clusters. 

More sophisticated measures of the morphology of X- 
ray contours are based on the moments of the squared- 
density distribution (Dutta 1995) or of the two-dimensional 
potential (Buote & Tsai 1995). These are more closely re- 
lated to the degree of substructure and so are more sensitive 
to the density parameter. Using the latter method, Buote & 
Xu (1997) claim that SCDM is inconsistent with ROSAT ob- 



Table 3. KS statistics for distinguishing power ratios within a 
0.75/i~^Mpc aperture: A is the maximum difference in cumu- 
lative fraction and Afo.os is the number of clusters required to 
distinguish the samples with 95 percent confidence using a 1- 
sidcd KS test. = 1-0 refers to the combined sample tCDM 
and SCDM. 



Samples 



P2 

A TVo.os 



^3 

A Nom 



Pa 
A ATq.os 



ACDM OCDM 
ACDM Ho = 1.0 
OCDM Ho = 1.0 



0.122 
0.206 
0.309 



125 
45 
20 



0.184 
0.115 
0.238 



55 
140 
35 



0.194 50 
0.113 145 
0.268 25 



servations of an, albeit incomplete, sample of bright clusters 
(Buote & Tsai, 1996). 

We define power ratios in a similar way to Buote & Tsai 
(1995): 

2 (a™ + 6^) 



Pm{Rc 



where 



R<Rc 



Rc 
R 



(6) 

(7) 

(8) 



Here (7?, 0) are projected polar co-ordinates about the clus- 
ter centroid , p is the density and the sum extends over all 
particles within 2 Abell radii. 

As might be expected the values of the power ratios are 
lowest for OCDM (least substructure) and greatest for the 
two critical density models, rCDM and SCDM. There is a 
slight tendency for high mass clusters to show more sub- 
structure, but this effect is smaller than the cosmological 
differences. The difference in the cumulative distributions of 
fm(0.75 Mpc) and the corresponding number of clusters 
required to distinguish the different cosmologies using a KS 
test are shown in Table ^; the critical-density models have 
been combined as they are indistinguishable. In contrast to 
Buote & Xu (1997), we find that P2 is best able to dis- 
tinguish low-density from critical-density models, whereas 
P3 and P4 do better at discriminating between ACDM and 
OCDM. Of course the numbers in Table ^ assume a direct 
correspondence between the theoretical distribution of dark 
matter and the observed X-ray flux. This association needs 
to be checked and a practical application of the method is 
likely to require more clusters than is indicated here. 

We have also calculated the first moment, Pi, within 
an aperture centred on the density maximum. This gives 
similar results to P2 because they both measure the dipole 
moment of the cluster. 

The different moments, or power ratios, are highly cor- 
related. The strongest relation, shown in Figure H, is be- 
tween P2 and P4. Interestingly the f2o = 1 cosmologies and 
OCDM show a tighter correlation than ACDM, that seems 
to extend to higher values of P4 for a given P2. This may 
be another way to discriminate between models with and 
without a cosmological constant. 

When we repeated the above tests for power ratios de- 
fined within an aperture of variable size, Rc = rigo, propor- 
tional to the virial radius of the cluster, then the differences 
between clusters in the ACDM and OCDM cosmologies were 
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Figure 9. P2 versus P4 in 0.75 h ^Mpc apertures for the different 
cosmologies. Only a subset of all clusters is sliown, for clarity. 



reduced, but the distinction between these and the critical- 
density universes was maintained. 

3.4 Density profiles 

Our final measure of the degree of isolation of a cluster is 
the slope of the density profile. We use the mean slope, s, 
between the radius, riso, within which the mean density is 
180 times the critical density (defined in Section i.3), and 0.1 
times this radius. Assuming a power-law density distribution 
p oc r"" in this range, then s is defined implicitly as 
^3-s _ f 

(9) 



m(0.32r) - m(O.lr) 10(3-'')/2 - 1 

where m(r) is the mass within radius r. 

s is well-defined only for clusters with a fairly smooth 
density distribution and so we restrict the analysis to those 
with Mint < 1.2. Histograms of the distributions of s for low- 
density and critical-density runs are shown in Figure ^[ 
The two are marginally distinguishable: assuming perfect 
data one would need a sample of 90 clusters to distinguish 
the two at 5 percent confidence. The mean slopes are 2.45 
and 2.40, respectively, consistent with our earlier conclusion 
that clusters in low-density universes are slightly more iso- 
lated (remember too that there are more clusters with sub- 
structure in the critical-density universes which have been 
omitted from the analysis). 

These differences are less than have been found in 
previous work. M95 found s to be much larger in open, 
low-density cosmologies than in flat ones, whereas Jing et 
al. (1995) found a weaker effect, but still larger than ours. 
We all agree that, for a given density parameter, clusters in 
open cosmologies are more isolated (i.e. have steeper density 
profiles) than those in fiat cosmologies. 

There is a weak trend of decreasing s with increasing 
mass in low-density cosmologies and increasing s with in- 
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Figure 10. Values of the slope of the density profile, as defined 
in Equation ^, for clusters with Mint < 1.2. The arrow shows the 
mean slope of the Navarro, Frenk & White (1995) profile in this 



creasing mass in critical-density cosmologies. Thus low-mass 
clusters show greater differences in their slope than do high- 
mass ones. As our simulations are of larger volumes, and 
hence contain more high-mass clusters than previous stud- 
ies, this may partly explain the smaller variation in s that 
we find. 

The mean values of s are similar to that of the density 
profile of Navarro, Frenk & White (1995), shown by the 
arrow in Figure hd. This will be discussed further in the 
next section. 

Using a fixed physical radius, rather than one defined 
in terms of overdensity is much less useful. Because in more 
massive clusters one is measuring the slope at a smaller frac- 
tion of the cluster radius, and because the density profile 
steepens with increasing radius, there is a decrease in s with 
increasing mass. The measured slope then becomes strongly 
dependent upon the mass of the cluster sample. 

It is inconceivable that the differences in slope, being 
much less than the spread in slopes even for the subset of 
smooth clusters, could be used to discriminate between dif- 
ferent cosmological models. 



3.5 Velocity dispersion profiles 

Figure ^ shows averaged velocity profiles for the clusters in 
cosmologies ACDM and SCDM (the others are similar). In 
each case we have included only relatively smooth clusters 
with Mint< 1.2 and have centred the profiles on the density 
maximum. Th e radii are scaled to the radius, rigo (defined 
in Section 3.3), at which the density drops to 180 times the 
critical density, and all clusters are given the same weighting 
in the average. For an isothermal profile ri8o is a fraction 
1/^3 « 0.58 of the virial radius (within which the mean 
overdensity is 180). 

The circles in the Figure show the normalised circular 
velocity profiles, Vc = (m riso/misor)"'^ where m is the mass 
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Figure 11. Velocity profiles for the combined sample of 100 
ACDM clusters with Mint< 1.2: circles — circular speed, solid 
line — fit to circular speed, dashed line — predicted circular speed 
from steady-state Jeans' Equation, stars — velocity dispersion, 
crosses — anisotropy parameter. 
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Figure 11 — continued Velocity profiles for the combined sample 
of 97 SCDM clusters with Mint< 1.2. Legend as for ACDM plot. 



contained within radius r and miso is the total mass within 
riso- These have been fit by a function of the form 



2 _ a;i8o(l + a;i8o)''"^ [(1 + xY ^ - (s - \)x - l] 
x (1 + a:)"-! [(1 + a;i8o)=-^ ^ - l)a^i80 - l] 
(s 7^ 2) corresponding to a density profile 
1 



p cc 



x{l + x)''' 



(10) 



(11) 



where x = r/a, xigo = '"iso/fl and the characteristic radius a 
and slope s are parameters of the fit. For s = 2, Equation |l^ 
must be modified slightly to 



Vc\s = 2) = 



a;i8o(l + a;i8o) [(1 + x) ln(l + x) — x^ 
x{l + x) [(1 + a;i8o) ln(l + a;i8o) - a;i8o] 



(12) 



Label 




{riso)/h ^kpc 


s 


a/r-180 


OCDM 


54 


0.62 


2.02 ±0.01 


0.29 ±0.04 


ACDM 


100 


0.56 


1.82 ± 0.01 


0.22 ±0.01 


rCDM 


65 


0.62 


2.01 ± 0.03 


0.37 ±0.01 


SCDM 


97 


0.57 


1.99 ± 0.02 


0.31 ±0.03 



s = 2 corresponds to the profile suggested by Navarro, Frenk 
& White (1995, hereafter NFW; 1996) to be a good fit to 
the galaxy clusters in their simulations, and s = 3 to the 
model proposed by Hernquist (1990) for the halos of ellip- 
tical galaxies. The low-density models were fit over a range 
0.1 < r/riso < 4 and the critical-density ones over a range 
0.1 < r/ri8o < 3. The lower radius in each case corresponds 
to about 60/i~^kpc, comfortably larger than the softening, 
and the outer radius was determined by the break in the ve- 
locity dispersion profile. The fits are shown in Table ^ and 
by the solid line in the Figure. 

In each case the extrapolated fit overestimates the cir- 
cular speed within r = O.lriso which implies that the inner 
density profile is shallower than 1/r; this may, however, be 
due to the finite softening which ranges up to 30 /i~^kpc for 
the critical-density runs. The clusters in all but one cosmol- 
ogy follow the NFW profile closely over the range in which 
they are fitted: the ACDM clusters would prefer a slightly 
shallower slope at r ~ ri8o but steepen at large radii to 
agree with the others. 

The 'concentration parameter', c « \/3ri8o/a, is 
around 6, in rough agreement with the values found by 
NFW for halos of this mass. The low-density cosmologies 
have higher values of c refiecting the higher redshift of clus- 
ter formation. 

The normalised velocity 

dispersion, a = (v'^riso/SGmiso)^^'^ , is shown by the stars 
in the Figure. In each case it rises gently out to r/riso ~0.2- 
0.3, then declines at larger radii. The variation in within 
the virial radius is about a factor of three; similar temper- 
ature variations are seen in hydrodynamic simulations (e.g. 
Thomas & Couchman 1992). Jing & Borner (1995) found 
a significant trend for steeper velocity dispersion profiles in 
open cosmologies when measured at constant physical radii. 
We suspect that this trend may disappear if the profiles are 
normalised to riso, as in Figure |l^. However, the effect is 
anyway quite weak for Qq = 0.3 and so we will not pursue 
it here. 

The variation of velocity dispersion with radius is not 
equal in each component. The anisotropy parameter, /3 — 
1 — CTt^/o"r^, where and at are the radial and (1-D) tan- 
gential components of the velocity dispersion respectively, 
is shown by the crosses in Figure ^ The orbits are ap- 
proximately isotropic in the core of the cluster but become 
increasingly radial as one moves out in radius. Contrary to 
the report in Crone, Evrard & Richstone (1994), we find no 
marked difference in the behaviour of the anisotropy param- 
eter between different cosmologies. 

It is usual to treat the central regions of clusters as in a 
quasi-steady state. We test this assumption by checking to 
see if the clusters satisfy the spherically-symmetric, steady- 
state Jeans Equation: 
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Figure 12. Values of the spin parameter for tlie combined sample 
of clusters. 



l_d(pcr,2) 
P 



dr 



+ 



(13) 



The value of Vc predicted by this equation is shown by the 
dashed line in Figure |ll| It is clear that the steady-state 
assumption is a reasonable one out to the virial radius. 



3.6 Angular momentum 

We measure the angular momentum of the clusters in terms 
of the dimensionless spin parameter, 



A 



where J is the magnitude of the angular momentum, T is the 
kinetic energy and m is the mass of the cluster. The original 
definition (Peebles 1971) uses the total energy of the system 
rather than the kinetic energy, but this is much harder to 
determine and the two are equivalent for a virialised system. 

The distributions of A are indistinguishable for each of 
our cosmologies and the combined sample is shown in Fig- 
ure |l|. The mean value of A is 0.060 (the 10, 50 and 90 
percentiles are 0.020, 0.051 and 0.108) but there is a sig- 
nificant trend of increasing mean and variance in A as one 
moves to smaller masses. 



4 CONCLUSIONS 

We have examined the structural differences between clus- 
ters in simulations of four cosmological models. Our basic 
conclusion is that relaxed clusters have very similar proper- 
ties in all cosmologies, but that there is a greater proportion 
of dynamically young clusters in critical-density universes 
than in low-density ones. It is hard to discriminate between 
low-density universes with and without a cosmological con- 
stant. 

The best way to constrain the density parameter is to 
look for substructure in the cluster population. Thus, in Sec- 
tion 3.2 it was shown that clusters in critical-density simula- 



tions are much more likely to have multiple cores, and that a 



very simple test, the centre-of-mass shift as the density con- 
tour is varied, can distinguish the two in a population of just 
20 clusters. The dipole moment of the squared-density (i.e. 
the P2 power ratio. Section 3.3) requires a similar sample 
size. 

Secondary measures of substructure, such as the shape 
of the iso-density contours, perform less well. The elliptici- 
ties we see (in agreement with other recent work) are sub- 
stantially smaller than M95 originally suggested. This may, 
in part, be due to our larger boxes. Most of the clusters 
which we find appear to be located at the intersections of 
large-scale filaments: in this sense none of them are truly iso- 
lated and many are still forming today even in low-density 
universes. 

The density and velocity-dispersion profiles of clusters 
have both been suggested in the pas t as cos molo gical probes. 
However, we have shown in Sections 3.4 and 3.5 that the pro- 
files of relaxed clusters are very similar in all models we con- 
sidered. They are in a quasi-steady state within the virial 
radius (that encloses a mean overdensity of 180), and are 
well-fit by the density profile of NFW. If we look at the com- 
plete sample of clusters, then the mean density and veloc- 
ity dispersion profiles decline less rapidly in critical-density 
models that are less isolated and have less well-defined cen- 
tres. However, we feel that this is a secondary effect and that 
it is better to look for substructure directly. 

Some of the properties that we have investigated (such 
as the degree of isolation, the major/minor-axis ratio and the 
angular momentum) have a slight mass-dependence. Others, 
(such as the density and velocity-dispersion profiles) show 
more regularity when measured in units of the virial radius, 
rather than a fixed physical size. This suggests that some of 
the variation in cluster properties which have been reported 
in earlier work may be due to a variation in the multiplic- 
ity function in different cosmologies. The redshift evolution 
of the number density of clusters is also expected to show 
a strong cosmological dependence. This two topics will be 
investigated in subsequent papers. 
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